New analysis of the SN 1987A neutrinos with a flexible spectral shape 
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We analyze the neutrino events from the supernova (SN) 1987A detected by the Kamiokande II 
(KII) and Irvine-Michigan-Brookhaven (1MB) experiments. For the time-integrated flux we assume 
a quasi-thermal spectrum of the form (_E/i5o)" e"'"^^'^''^'' where a plays the role of a spectral 
index. This simple representation not only allows one to fit the total energy Etot emitted in Ve 
and the average energy {Ey^), but also accommodates a wide range of shapes, notably anti-pinched 
spectra that are broader than a thermal distribution. We find that the pile-up of low-energy events 
near threshold in KII forces the best-fit value for a to the lowest value of any assumed prior range. 
This applies to the KII events alone as well as to a common analysis of the two data sets. The 
preference of the data for an "unphysical" spectral shape implies that one can extract meaningful 
values for {Ej^^) and -Etot only if one fixes a prior value for a. The tension between the KII and 
1MB data sets and theoretical expectations for {Ei^^ ) is not resolved by an anti-pinched spectrum. 

PACS numbers: 14.60.Pq, 95.55.Vj, 95.85.Ry, 97.60.Bw 



I. INTRODUCTION 

The neutrino observations of supernova (SN) 1987A 
[1-6] have long been taken as a confirmation of the 
salient features of our physical understanding of the core- 
collapse SN phenomenon. At the same time, a detailed 
interpretation [7-14] of the data is difficult because of 
a number of "anomalies" [15]. In particular, the energy 
spectra implied by the Kamiokande II (KII) and Irvine- 
Michigan-Brookhaven (1MB) events are barely consistent 
with each other. Moreover, the average v,, energy implied 
by KII is much lower than expected from numerical sim- 
ulations. The energies implied by 1MB alone would be 
compatible with theoretical models, but the results of 
a common KII and 1MB analysis are difficult to square 
with expectations. 

In the presence of neutrino oscillations, the v,, flux ob- 
served in a detector is a superposition of the v,, and i>^.r 
fluxes produced at the source. Therefore, the KII and 
1MB detectors may have observed different fluxes be- 
cause the SN I987A neutrinos have traversed different 
sections of the Earth so that matter effects can mod- 
ify the oberved superposition of source spectra [16, 17]. 
However, it is expected that (Sp^^) > {E^j^) so that 
oscillations aggravate the tension between observed and 
expected energies [11, 16, 17]. In any event, the dif- 
ferences between (iJp^ ^ ) and (Sp^ ) are probably much 
smaller than had been thought previously [28-30] so 
that the Earth matter effect is no longer expected to 
cause gross modifications of the observed Ve spectrum. 
Of course, the relatively subtle modifications caused by 
Earth matter effects can be crucial for identifying the 
neutrino mass ordering from the high-statistics neutrino 
signal of a future galactic SN [18-27]. 

Previous studies of the SN 1987A neutrinos usually as- 
sumed a thermal spectrum and then extracted ) and 
the overall flux from the individual detector signals or 
from a combined analysis. One exception is the analysis 



of Janka and Hillebrandt [9] who assumed an effective 
Fermi-Dirac distribution of the form 



Lp{E) oc 
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where T is an effective temperature and 77 a degeneracy 
parameter. In their maximum-likelihood analysis they 
allowed only for positive values of 77. With this prior 
they found a best-fit value of 77 = for both the KII and 
1MB data sets, suggesting that the data prefer the broad- 
est allowed distribution compatible with the prior range 
for 77. However, even allowing for negative values for 77 
would not have changed these results because 77 ^ — 00 
corresponds to a Maxwell-Boltzmann spectrum, differing 
only marginally from the Fermi-Dirac case with 7; = 0. 

The main purpose of our new study is to investigate if 
better internal agreement of the SN 1987A data as well as 
better agreement between the data and theoretical expec- 
tations can be achieved if a more flexible representation of 
the spectral shape is assumed. Numerical studies of neu- 
trino transport suggest that the instantaneous spectra 
are "pinched," i.e. that they are narrower than a thermal 
spectrum [9, 29]. A pinched spectrum can be represented 
by a Fermi-Dirac distribution with positive 77. However, 
the SN 1987A data measure the time-integrated spec- 
trum, i.e. a superposition of instantaneous spectra with 
varying average energies. Therefore, the integrated spec- 
trum may well be broader, not narrower, than a thermal 
spectrum, i.e. it may well be anti-pinched. 

If the time-intergrated spectrum is quasi-thermal in 
the sense that it rises from zero for low energies, reaches 
a maximum, and has a long tail to high energies, then 
the simplest conceivable representation is [29] 
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1 and Eq is an energy scale with the 
The numerical parameter a controls 
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We note that a = 2 corresponds to a Maxwell-Boltzmann 
spectrum, a > 2 to a pinched spectrum with suppressed 
high- and low-energy tails, and a < 2 to an anti-pinched 
spectrum. For a — > oo the spectrum becomes 5{E — Eq). 
Since the oscillations effect on the detected spectrum is 
reasonably small (sec Sec. I), it makes sense to fit directly 
the effective Ue spectrum, after the oscillations, with the 
distribution ofEq. (2). 

In Sec. II we present the SN 1987A neutrino data and 
perform a new maximum-likelihood analysis. A summary 
and conclusions are given in Sec. III. 



II. MAXIMUM-LIKELIHOOD ANALYSIS 

A. SN 1987A Data 

We limit our analysis to the SN 1987A data of the 
KII [1, 2] and 1MB [3, 4] detectors that measure SN 
neutrinos almost exclusively by the inverse beta reaction 
De+P ^ n + e'^ . We show the measured positron spectra 
in Fig. 1 where we have left out the KII event No. 6 that 
is attributed to background. In Table I we summarize the 
properties of the neutrino signal in the two detectors. We 
do not include the Baksan Scintillator Telescope (BST) 
data [5, 6] because it is much more uncertain which of 
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FIG. 1: Positron spectra detected at KII (upper panel) and 
1MB (lower panel) in connection with SN 1987A. We also 
show the energies of the individual events together with their 
experimental errors. 



TABLE I: Number of SN 1987A events and average positron 
energies at the KII and 1MB detectors. 



Detector 


events 


) [MeV] 


KII 


11 


15.4 ± 1.1 


1MB 


8 


31.9 ±2.3 



the events have to be attributed to background, i.e. in 
a maximum-likelihood analysis one would have to model 
the background. This requires to use the time structure 
of the neutrino burst [13, 14], whereas we limit our study 
to the time-integrated flux. 



B. McLxwell-Boltzmann Spectrum 

We perform a maximum-likelihood analysis of the 
SN 1987A signal along the lines of the previous litera- 
ture such as Rcf. [11]. For the detection cross section 
we use an updated analytic fit [31] and for the detection 
efficiencies we use the analytic fit functions of Ref. [32]. 

In order to compare with the previous literature we 
first consider a MaxwcU-Boltzmann i>e spcctriim, i.e. we 
use a = 2 in Eq. (2). As fit parameters wc use the av- 
erage De energy Eq and the total energy E^tot emitted by 
SN 1987A in the form of Pg- Of course, these parameters 
refer to the spectrum measured in the detectors after the 
partial flavor swapping caused by neutrino oscillations. 
Our results shown in Fig. 2 agree with the previous lit- 
erature and illustrate once more the tension between the 
average energies implied by the two detectors. 
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FIG. 2: Contours of constant likelihood that correspond to 

68.3%, 95% and 99% C.L. in the plane Eo-Etot for an assumed 
Maxwell-Boltzmann i^e spectrum, corresponding to a = 2 in 
Eq. (2). The dotted hnos refer to the KII and 1MB signals, 
respectively, whereas the solid lines represent a joint analyis. 
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C. General Spectrum 

Next we use a general spectrum of the form Eq. (2) 
with Eq, Etot, and a as our fit parameters with a prior 
a > 0. We show the best-fit values for these parameters 
as well as the implied event numbers and average positron 
energies in the detectors in Table II. This analysis is 
performed for each detector separately as well as for the 
joint case. 



TABLE II: Best-fit values for Etot in lO''^ erg, ^Eo in MeV, 
and a based on the indicated data sets. The implied charac- 
teristics of the expected detector signals are also shown. 





Best-fit param's 


KII 




1MB 




Data Set 


-Etot Eo a 


events 




events 




KII 


17.4 3.7 0.0 


10.9 


15.3 


1.3 


27.5 


1MB 


1.1 30.6 60.0 


6.6 


30.3 


7.9 


31.4 


Joint 


10.1 5.4 0.0 


14.1 


19.2 


4.9 


32.0 



The 1MB data alone prefer a functional form with a 
very large value for a, i.e. essentially a S function. This 
behavior is intuitively obvious because 5 of the 8 1MB 
events have positron energies in the very narrow range 
36-39 MeV. The KII data alone, on the other hand, pre- 
fer a = 0, i.e. a huge total flux of low-energy Dg. The pile- 
up of events just above the KII threshold of 7.5 MeV im- 
plies that the turn-over at low energies of a quasi-thermal 
spectrum is not visible in the data, which are sensitive 
only to the decreasing higher energy tail of the spectrum. 
This feature is reproduced in the likelihood by pushing 
Eq toward low values and broadening as much as possible 
the width. Moreover, the low value of Eq is compensated 
by a huge amount of the total emitted energy Etot ■ We 
have checked that even when we allow for negative a val- 
ues the maximum of the likelihood always coincides with 
the lowest allowed value. Of course, for a < —1 the total 
neutrino flux diverges. 

Contrary to our expectation, this behavior remains the 
same for a joint analysis of the two data sets, even though 
the total flux is smaller and the average Ve energy is 
larger. 

In Fig. 3 we show the difference A In £ relative to the 
best-fit value as a function of a for the two separate data 
sets and the joint analysis. In all cases we have marginal- 
ized over the parameters Eq and i^tot- The 95% C.L. 
allowed range for a single degree of freedom corresponds 
to Aln£ = 1.92 shown as a horizontal dot-dashed line 
in Fig. 3. Evidently 1MB alone has no strong preference 
for any spectral shape. The 95% upper limit on a is 
lower for the combined data than for KII alone because 
the combined data also prefer a larger Eq and lower to- 
tal flux. Overall, the data do not distinguish in a useful 
way between different plausible spectral shapes. In par- 
ticular, the data do not prefer a quasi-thermal spectrum 
but rather a monotonically falling one, in contrast with 
plausible expectations. 
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FIG. 3: Curves of Aln£ relative to the best-fit value as a 
function of a for KII (dotted fine), 1MB (dashed line) and 
for the joint data set (continuous line). The 95% C.L. for 
a is indicated with the dot-dashed horizontal line. We have 
marginalized over Eq and -Etot- 



D. Fixed Prior Values for a 

We conclude that in order to extract information on 
the average Pe energy and the total flux it is not useful 
to marginalize over the parameter a because the result 
depends sensitively on the assumed prior range for a and 
because unphysically low values are prefered by the data. 
Therefore, we return to fixing the spectral shape in ad- 
vance and analyze the data for different choices of a. In 
Fig. 4 we show the 95% C.L. contours in the EQ-Etot 
plane for a = 0, 2 and 4 where the case a = 2 is identical 
with the 95% C.L. shown in Fig. 2. Larger values of a 
(pinched spectra) increase the range of average energies 
and bring them closer to theoretical expectations, but at 
the same time decrease the range of overlap between 1MB 



TABLE III: Best-fit values for -Etot in 10^^ erg and Eo in MeV 
for the indicated fixed choices of a. The implied characteris- 
tics of the expected detector signals are also shown. 



Best-fit KII 1MB 







Data Set 


Etot 


Eo 


events 




events 




a 


= 


KII 


17.4 


3.7 


10.9 


15.3 


1.3 


27.5 






1MB 


23.7 


5.0 


28.7 


18.3 


8.0 


30.8 






Joint 


10.1 


5.4 


14.1 


19.2 


4.9 


32.0 


a 


= 2 


KII 


8.9 


7.9 


11.0 


15.3 


1.1 


26.7 






1MB 


8.2 


11.6 


20.6 


20.1 


8.0 


30.9 






Joint 


5.9 


11.2 


14.0 


19.5 


4.9 


30.3 


a 


= 4 


KII 


6.6 


10.2 


11.0 


15.4 


1.0 


25.3 






1MB 


4.7 


15.9 


16.3 


21.8 


8.0 


31.0 






Joint 


4.7 


14.2 


14.1 


19.8 


4.8 


29.5 



4 
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FIG. 4: Contours of 95 % C.L. in the plane of Eo and iJtot 
for a = (top panel), a = 2 (middle) and a = 4 (bottom). 
The dotted lines are for KII and 1MB separately whereas the 
solid line is for the joint data set. 



and KII. In Tabic III wc show the best-fit parameters Eq 
and i^tot for the three a cases and different data sets as 
well as the implied signal characteristics. 



III. DISCUSSION AND SUMMARY 

The expected neutrino energies from core-collapse SNe 
are of practical importance for sensitivity forecasts of a 
future galactic SN neutrino detection, in particular with 
regard to neutrino mixing parameters [18 27]. The ex- 
pected energies are also important in the context of cur- 
rent limits and future detection possibilities of the diffuse 
SN neutrino background from all SNe in the universe [33- 
37] . The internal tension within the SN 1987A neutrino 
data as well as the tension with theoretical expectations 
has persuaded most workers in this field to ignore the 
data and rely on the output of numerical simulations 
even though current SN theory may still be missing an 
important piece of input physics to obtain robust explo- 
sions [38]. 

One explanation for the tension between the KII and 
1MB data is that the detectors actually observed different 
spectra due to the Earth matter effect [16, 17]. However, 
our current understanding of flavor-dependent neutrino 
spectra formation suggests that the flavor-dependent av- 
erage energies in the anti-neutrino sector are not very 
different [28, 29]. Moreover, if the observed 9e had been 
born as higher-energy 9^,^ at the source, the tension be- 
tween theoretically expected and actually observed i>e en- 
ergies would be worse. 

Loredo and Lamb [14] have tested a large variety of 
neutrino emission models and in particular have included 
what amounts to a bi-modal spectral shape that con- 
tains a lower-energy component attributed to the SN ac- 
cretion phase and a higher-energy one attributed to the 
neutron-star cooling phase. They stress that such a spec- 
tral form is strongly favored by the data relative to a 
single-mode thermal spectrum. However, numerical sim- 
ulations do not predict a bimodal spectrum because the 
average energies continuously increase from the accretion 
to the neutron-star cooling phase and at late times de- 
crease again. Therefore, the spectral shape of the time- 
integrated flux does not seem to exhibit a bi-modal form 
but rather is expected to be a broadened quasi-thermal 
spectriim. 

Motivated by this observation we have analysed the 
SN 1987A data, assuming a quasi-thermal spectrum of 
the general form Eq. (2) that is flexible enough to ac- 
commodate a continuum of spectral shapes from narrow 
(pinched) to broadened (anti-pinchcd) spectra relative to 
a Maxwell-Boltzmann distribution. 

Perhaps unsurprisingly in view of the tension between 
the KII and 1MB data, wc find that the broadest possi- 
ble distribution allowed by the chosen prior range of a is 
prefered. In this way the tension between the data sets 
is somewhat reduced, but at the same time the average 
i/e energies are pulled to lower values, thus exacerbat- 
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ing the tension with the output of representative numer- 
ical simulations. (The average neutrino energies found in 
many different numerical simulations have been collected 
in Ref. [29].) 

Assuming the time-intergrated neutrino flux from a SN 
indeed exhibits a quasi-thermal spectrum roughly of the 
form Eq. (2), the implied average t'e energies and total 
emitted energy depend sensitively on the chosen prior 
range for a. The data themselves prefer the smallest 
possible a-values, i.e. not a quasi-thermal spectrum but 
rather a monotonically falling one. The assumption of 
a realistic quasi-thermal spectrum is not simultaneously 
consistent with typical theoretical expectations for {Ep^ ) 
as well as the separate data sets from 1MB and KII. 
Therefore, it remains unresolved if the SN 1987A data 
or theoretical simulations give us better benchmarks for 
the average Pe energies to be used, for example, in the 
context of searches for the cosmic diffuse SN neutrino 
background. 



It appears that the question of the true neutrino spec- 
trum from a typical SN can be empirically resolved only 
by the high-statistics signal from a future galactic SN or 
by the patient accumulation of data on a neutrino-by- 
neutrino basis from SNe in nearby galaxies [39]. 
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